Transitions to equilibrium state in classical 4 lattice 
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Abstract 

> 

I Statistical behavior of a classical (/> 4 Hamiltonian lattice is investigated from 

o . 

microscopic dynamics. The largest Lyapunov exponent and entropies are 

o 

considered for manifesting chaos and equipartition behaviors of the system. 
It is found, for the first time, that for any large while finite system size there 
exist two critical couplings for the transitions to equipartitions, and the scaling 
behaviors of these lower and upper critical couplings vs the system size is 



> 

X 
b . 

numerically obtained. 
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I. INTRODUCTION 



The investigation of statistical behavior of Hamiltonian coupled oscillators from micro- 
scopic dynamics has attracted much attention since the original Fermi-Pasta-Ulam (FPU) 
problem [1-4]. The important issue in this aspect is the approach to equilibrium of non- 
linear system with many degrees of freedom. It was found in the FPU study that with 
sufficiently small energy level various modes of the system do not go to equipartition state 
as conventionally expected for equilibrium [1]. This observation has stimulated many works 
studying the transitions from nonequilibrium to equilibrium, based on numerical simulations 
of microscopic dynamics [5-9]. Various models of Hamiltonian system with many degrees 
of freedom have been considered for this purpose, among which the 4 lattice is one of the 
most extensively investigated . 



The model Hamiltonian reads 



i=i,....N 2 2 4 

Qn+i = <h , Pn+1 = Pi 



With the rescalings % = y A*/ Pi = y m /f3pi and t = Jm/nt, we can reduce Eq. (1) to 

2 / \2 4 

H — i^o + t + ~ >r ~k)- So far most of papers considering the statistical behavior of 

- I Y 

the lattice 4 model have used this form and taken the average initial energy as the control 
parameter. In this paper, we use an alternative while equivalent method to fix the initial 
average energy E = ■h Ei and vary the coupling strength fi as the control parameter, 

i=l,. ..N 

i.e., we consider the lattice 4 Hamiltonian 

H= £ 4 + + (2) 

i=l,....N Z Z 4 

which leads to the set of canonical equations 

dH 3 
Pi = —a— = ~Qi + A*(?<+1 + ft-l - 2 <?i) 
Oft 

?< = a" = ft 3 



We take E — 1.25 x 1CT 3 throughout the paper without lossing any generality. 

Before going to the detail of the dynamics of Eqs. (3), we can first briefly discuss the 
general feature of the system for different coupling regimes. For /i — > 0, the system consists 
of practically independent oscillators, each has one degree of freedom, and thus has nearly 
periodic motion. The dynamics of the system may be quasiperiodic and the approach to 
equilibrium is not possible. For fi — > oo, the nonlinear terms qf/4 become not important in 
comparison with the coupling. Thus, different modes are practically uncoupled to each other 
[10-12], and the system dynamics is quasiperiodic again. Complicated dynamics approaching 
equilibrium may be expected for not too small and not too large couplings (i.e., not too 
large and not too small energies if we fix ji = 1 and use the average energy as our control 
parameter). Up to date, all papers discussing the statistical behaviors of the 4 model 
consider only small energy (i.e., large coupling /x) regime. To our knowledge, the behaviors 
of large energy (or say, small coupling jj) regime has never been discussed. In this paper 
we will study the microscopic dynamics of the system and the related statistical features in 
the whole coupling range by paying particular attention to the new area, the small coupling 
regime. 

In the following, we will make a Lyapunov exponent analysis of the system dynamics in 
Sec. II; and study the statistical quantities, the entropy and the probability distribution of 
system variables in Sec. III. In Sec. IV, we will consider the scaling behaviors observed in 
the transitions to equilibrium, and some conclusions and possible future developments will 
be discussed in the last section. 

II. LYAPUNOV EXPONENT ANALYSIS 

Given the average energy E and the coupling /i, the dynamics of Eqs. (3) is determined 
by the initial conditions of the (%,£>«) variables. In this paper we will test and compare two 
types of typical initial conditions. First, the initial energy is distributed to a single k = 1 
mode 
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Pi = Asm(2iri/N) + er» , \(Ji\ < 1 , A 




(4) 



z = 



1,2, ...,N 



i.e., 



, k = 0,1, 2,. ..JV/2 



Pi 



A//V/2, P fc = , k = 0, 2, 3, N/2 



Qk(t) 



2/N X! cos(2Trki/N), P k (t) =Q k (t) 



i=i 



This initial condition has been popularly used in studying the statistical behavior of Hamil- 



which is considered for the first time, to our knowledge, in this paper. Here the equipartition 
problem will be analysed for both mode distribution and site distribution. For avoiding the 
strict symmetry restriction preventing equipartitions we add a small random numbers G\ in 
Eqs. (4) and (5) with \a,\ < 1(T 7 . 

Lyapunov exponent (LE) analysis is a very useful dynamical tool for revealing the sta- 
tistical behavior of Hamiltonian systems [10,12,13]. If the dynamics is quasiperiodic, the 
largest LE is zero, and the system will never reach the equal-probability-distribution in the 
energy surface for N > 2. Because for a Hamiltonian system with N degrees of freedom the 
quasiperiodic trajectory can move only in a iV-dimensional torus, while the energy surface 
has 2N — 1 dimensions. Therefore, the positivity of the largest LE is a strong indication (of 
course, not a sufficient condition) of possible equilibrium state. 

It is well known that there exist some difficulties for correctly computing LEs in Eqs. 
(3) when iV is large [10,12]. Numerically, LE is calculated approximately with a very long 
but still finite time T. However, if iV is large, it may happen that the system wonders in a 



tonian coupled oscillators flOfl . Second, the initial excitation is restricted on a single site 



qi = 0, z = l,2,....,JV 




(5) 



certain A^-dimensional torus for an extremely long time and then bursts to a chaotic status. 
In Figs. 1(a) and (b) we take N = 64 and use the k — 1 mode initial condition (4), and plot 
the largest LE Ai against the average time T for fi = 6.9 x 10~ 7 and 2.5, respectively. In 
both case we find that Ai is maintained at zero for long time, and jumps to be positive at 
a certain large T. If we decrease \x from \x = 6.85 x 10~ 7 and increase \x from \x = 2.9, the 
periods for Ai to change from zero to positive increase. It seems that there exist two critical 
couplings Hi and /12, and zero Ai is observed for \i < ^i or // > and chaos is identified 
between these two couplings fii < \i < Hi- Transitions from quansiperiodicity to chaos occur 
at both large (^2) and small (/ii) critical couplings. For approximately fixing \i\ an d ^2, we 
take average time T = 1 x 10 7 , and zero Ai is assumed if no burst to positive LE accurs 
within the period T . 

In Fig. 1(c) we plot Ai vs /x of Eqs. (3) at = 32 and for both initial conditions 
(4) (circles) and (5) (triangles). In each curve we find indeed critical couplings and 
/x 2 (/4); an d Ai is zero for < fii (// < Hi) an d A* > A*2 (a* > ^2)- Positive Ai, indicating 
chaos, is observed between these two critical couplings. Another interesting point is that 
the circle and triangle lines have practical identical positive Ai in the region \i x < ji < 0.25, 
this verifies approach to equilibrium, which is independent of the initial conditions. 

In Fig. 1(d) we do the same as (c) by taking larger system N = 64. For the initial 
condition (4), the lower (higher) critical coupling Hii^) becomes much smaller (much larger) 
than that in (c). This feature is intuitively expected since it is familiar to us that for large 
one needs only small interactions between the subsystems to bring the system to equilibrium 
states fL2] . [note, for Eqs. (3) large coupling // corresponds to small interactions between 
different modes]. However, for the initial condition (5), the critical coupling h'i increases as 
the system size increases. This phenomenon, which seems to be against the intuition, can be 



easily understood, based on the measue synchronization (MS) of Hamiltonian systems |T4 
In order to reach the equilibrium state, all the oscillators of Eqs. (3) should have identical 
statistical behavior, independing on their initial condition, because they are identical in 
microscopic dynamics. The necessary condition realizing this identity is that the trajectories 



of all oscillators should be identical in measure, so-called MS. The critical coupling for the 
MS between any two neighbor oscillators depends on their energy difference. Large coupling 
is needed for the MS if a pair of neighbor oscillators have large energy difference. Since 
for condition (5) the energy difference between the site i — 1 and its neighbors increases 
by enlarging the system size with fixed average energy [in (5) E\ — E 2 = NE], the critcal 
coupling for MS, which is nothing but increases as well for increasing the system size 
(actually n' t oc y/~N). This feature is essentially different from the system behavior with the 
initial condition (4). 



III. ENTROPIES AND PROBABILITY DISTRIBUTIONS 

The dynamic behavior of LE is closely related to the statistical behavior of the system. 
It is interesting to further investigate the energy distribution of subsystems, because the 
equipartition behaviors of energies of different modes and different sites are some necessary 
conditions of the equilibrium state. 

The energies for different modes can be computed by 

E m (k,t) = ^(t) + ^lQl(t), fc = 0,2,....y (6) 

where uik = 2sin(^). At the same time, the energy for different sites are represented as 

E s (i,t) = ^ 2 (t) + \qf(t), i = 1,2, N (7) 

For the long time behavior, we are interested in the averages 

1 f T 



E m (k,T) = -( E m (k,t)dt 
1 Jo 



E a {i,T) = ^ F E.{i,t)dt (8) 
1 Jo 

In actual simulations we take large while finite T = 10 7 . 

Based on the definition of Eqs. (8), we can investigate the system statistics by defining 
entropies [IS| in the mode space 

6 



and in the site space 



N/2 

S m = -Y,n k \nQ k (9) 

Ai=0 

N/2 

n k = E m (k)/ E E m (k') 

k'=0 



N 

s s = -J2ti l \nn' t (io) 
i=i 

n; = ^(<)/ E ^.(0 
»'=i 

In numerical computation we measure the scaled entropy differences 

Vm = (S m (msx) - S m )/S m (max) 

7] s = (5s (max) - S s )/S s (max) (11) 

with S m (max) = ln(iV/2 + 1), ^(max) = ln(iV). In Fig. 2(a) [(b)] we plot i] m (r) s ) vs \i from 
the initial condition (4) [the initial condition (5)]. We find that in (a) i] m has discontinuous 
jumps at the two critical couplings /ii and /i 2 and it takes very small value for /ii < /i < /i 2 , 
indicating large S m entropy in this coupling regime. In (b) i] s has a plateau of i] s ss 1 for 
/i < and it jumps down to about r] s < 1 x 1CT 6 , and keeps this small value in the entire 
regime // > jj! v It is noticed that the second transition at fi' 2 is not manifested by rj s , because 
the MS between sites is not affected by this transition. 

In order to show that in the regime // : < /i < 1.65, the system actually reach the 
statistical equilibrium state for both initial conditions, we take /i — 0.1 and plot in Figs. 
2(c) and (d) the probability distributions of x and p variables, respectively, for all the 64 
particles and for both initial conditions (4) and (5). Each probability distribution curve 
is drawn with 1 x 10 6 data. From the figures it is verified that all the sites have identical 
distributions, showing perfect equipartion behaviors. 
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IV. SCALING BEHAVIORS IN THE TRANSITIONS TO EQUILIBRIUM 



In the above numerical analysis we found two transitions of Eqs. (3) for the initial condi- 
tion (4) from quansiperiodicity to chaos, one at small coupling transition to equipartion 
between different sites in real space, the other at large coupling /i 2 , transition to equiparti- 
tion between various modes in Fourier space. Both \i\ and /i 2 can be considerably changed 
by varying the system size. It is interesting to investigate the possible scaling behavior 
between the critical /1 and the system size N. 

In Fig. 3(a) we plot both ln/Z! and ln/i 2 vs IniV. Between the two curves the system 
motion is chaotic, where statistical equipartition behavior is expected. And quasiperiodicity 
exists in the regions below the fii curve and above the /i 2 curve. By quasiperiodicity we 
mean that quasiperiodic motion is observed in the time period t < T = 1 x 10 7 . And we 
assume, (without proof) that the actual critical couplings for the transitions are very near 
the plotted curves of Fig. 3(a). 

Both /ii(iV) and ^{N) curves have well behaved power law scalings 

a ^3.82 (12) 

/i 2 (iV) oc p « 1.92 (13) 

The solid lines fii(N) and fi 2 {N) are given by A/N a , A « 1 x 10 -3 , a « 3.82 and /J, 2 (N) = 
BN?, B « 2.4, (3 « 1.92, respectively. Equations (12) and (13) are numerical results. We 
do not know how to derive the exponents a and (3 analytically, and even have no idea at 
present to explain them. We suggest these scaling may be model dependent and do not have 
universal behavior. For instance, the scaling (13) is different from the scaling Casetti and 
coworkers obtained for the FPU model [I2|. 



With the initial condition (5) the critical point [i x and /i 2 are also related to N. The 
scalings are shown in Fig. 3(b), which read 

n[(N) oc iV 7 , 7 w 0.57 (14) 

/4(iV) oc N? , » 1.95 (15) 



V. CONCLUSION AND DISCUSSION 



In conclusion, we have investigated the statistical behavior of the classical 4 lattice 
from its microscopic dynamics. The particular attention focuses on the transitions to chaos, 
equipartition behavior and equilibrium state. By varying the coupling strength, two tran- 
sitions to equipartition from both small and large couplings are observed. The scaling 
behaviors of these couplings against the system size are revealed numerically. 

However, there are still some problems unsolved. First, it is difficult to precisely fix the 
transition points, because near the transitions the system may show zero LE for extremely 
long period, and suddenly jump to positive LE after a certain time. Thus, by a finite time 
simulation, one cannot guarantee zero LE, and consequently cannot definitely determine the 
critical conditions. In this paper, we apply an time truncation in identifying fii and /i 2 , i.e., 
we regard the largest LE to be zero if its zero value is kept for a given very long while still 



finite time [see Ref. JT2[]- We guess that the actual transition points are not far from the 
critical couplings shown in Figs. 1-3. We will try a systematical method to convincingly 
identify the transitions in the future work. Moreover, the mechanism underlying the scaling 
relations (12), (13) and (15) are still unclear, and worthwhile further investigating. 
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Captions of figures 

— N 

Fig. 1 The largest LE of Eqs. (3) Ai for the fixed average energy E = jr E% — 

i=i 

1.25 x 10~ 3 . (a) The largest LE vs the time T with T being the time for computing Ai. 
pi — 6.9 x 10~ 7 , N = 64, and the initial condition (4) is taken. It is shown that Ai(T) 
remains zero for a long period, and it bursts to be positive after a certain time, (b) the 
same as (a) with /j, = 2.5. (c) Ai vs /i for N — 32. Circles are for the initial condition (4) 
and triangles are for the initial condition (5). The transition points /ii, /j^ and fj! 2 are 
determined by computations of T = 1 x 10 7 . (d) The same as (c) with N = 64. 

Fig. 2 N = 64. i] m [(a)] and i] s [(b)] defined in Eq. (11) plotted vs ji. In (a) the initial 
condition (4) is used and in (b) Eq. (5) is used. i] m has discontinuities at two transition 
points Hi and /i 2 of Fig. 1(d), and r] s jumps down at the transition point \i x of Fig. 1(d). 
(c), (d) fj, = 0.1. The [(c)] and Pi [(d)] probability distributions of the sites i — 1,2, ....N 
for both initial conditions (4) and (5). 

Fig. 3 (a) The transition couplings /ii and /jl 2 plotted for different systems size with the 
initial condition (4) being taken. In the region between pii and \i2 lines, the systems motion 
is chaotic, and equipartition behavior can be identified, /ii and pi2 have power law relations 
with the size iV as /ii oc 1/iV 3 - 82 , /j, 2 oc A^ 1 - 92 . The exponents are measured by the slopes of 
the two solid lines, (b) The same as (a) with the initial condition (5) used. pi x oc A^ 57 , 
/4 oc N 1 - 95 . 
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